Carbon antisite clusters in SiC: a possible pathway to the Dn center 
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The photoluminescence center Dn is a persistent intrinsic defect which is common in all SiC 
polytypes. Its fingerprints are the characteristic phonon replicas in luminescence spectra. We 
perform ab-initio calculations of vibrational spectra for various defect complexes and find that carbon 
antisite clusters exhibit vibrational modes in the frequency range of the Dn spectrum. The clusters 
possess very high binding energies which guarantee their thermal stability — a known feature of the 
Dn center. The dicarbon antisite (C2)si (two carbon atoms sharing a silicon site) is an important 
building block of these clusters. 

PACS numbers: 61.72.-y, 63.20.Pw, 78.55.-m 



C/2 



X3 
O 

> 



0^ 
O 
cn 

o 



I 

O 

o 



I. INTRODUCTION 

Many unique features of SiC, such as the wide band 
gap and the very high electrical and thermal stability, 
render this semiconductor especially important for high 
power, high frequency and high temperature applica- 
tions. As for any semiconductor material, the identifi- 
cation and understanding of structural defects and im- 
purities is the key to the technological control of SiC. 
This is especially important in connection with the ion- 
implantation, which is inevitably accompanied by the 
generation of intrinsic defects. Some of these defects 
are thermally very stable. Two important examples are 
the photoluminescence (PL) centers Dj (Ref. Q) and Dn 
(Ref. 13). Although the centers have been extensively 
studied over the past 30 yearSfi*2i2i^iSiSii their microscopic 
origin remained unclear. In this work we study carbon 
antisite clusters and their possible relevance to Dn-type 
centers. 

Dii centers are always present in ion-implanted and an- 
nealed SiC samples, regardless of the implanted species 
and the polytype. They were also found (in low con- 
centration) in as-grown material.- The abundance of Dn 
centers increases significantly at high annealing temper- 
atures (above 1300°C) and the center is thermally stable 
up to 1700°C.'' These properties imply that the Dn cen- 
ters are intrinsic defect complexes. The Dn luminescence 
originates from the recombination of an exciton bound 
to the defect. The spectra are usually identified by the 
main zero-phonon line. The important specific feature 
of Dn are characteristic phonon replicas above the SiC 
phonon spectrum. In the original experiments in 3C- 
SiC the strongest five replicas were highlightedi^ In later 
experimentsSii in 4H and 6II-SiC additional localized vi- 
brational modes (LVMs) were observed (more than 12 
LVMs have been counted in 4H-SiG&). Many of these 
observed modes are polytype- independent. A compari- 
son of the spectra from Ref. 7. and the earlier Ref. jj re- 
veals that more than the five highlighted LVMs may also 
be contained in the 3C spectrum. Their spectral density 
resembles the phonon density of states in diamond, indi- 
cating that a carbon-dominated defect is responsible for 
the observed spectrum. 




FIG. 1: Structure of the dicarbon antisite (C2)si and a carbon 
cluster with four atoms on a silicon site (C4)si in 3C-SiC. 



In this work we consider a number of carbon-related 
defects in SiC: the carbon di-interstitial, which was the 
first suggested model for the Dn center )S the carbon anti- 
site, the carbon split-interstitials and several carbon clus- 
ters. Comparing the calculated LVMs with the experi- 
mental Dn spectra we are able to unequivocally rule out 
all candidates except the carbon antisite clusters. The 
latter are the only defects that provide a rich LVM spec- 
trum similar to that of the Dn center. The core struc- 
ture of all these clusters is the dicarbon antisite (C2)si 
(cf. Fig. n left) which acts as an aggegation center for 
further carbon atoms. In this paper we will first present 
our results for the dicarbon antisite and then consider 
the formation of more complex defect clusters. We calcu- 
late the LVMs associated with these defects and compare 
them with the Dn spectrum. Although theoretically all 
these defects can acquire different charge states, we fo- 
cus on the neutral state as the most relevant for PL. For 
excitons bound to a charged defect more non-radiative 
recombination channels are open, and they are unlikely 
to be seen in PL-experiments.^ We also find that the for- 
mation kinetics of the carbon clusters agrees with the 
known features of the Dn center. These results indicate 
that although neither of the considered defects can alone 
explain all features of the Dn center, the carbon antisite 
clusters may serve as a core structure of Dn-type defects. 



II. METHOD 

We employ an ab initio density functional theory 
(DFT) approach as implemented in the software package 
FHI96SPIN J. Smooth norm-conserving pseudopotentials 
of the TrouUier-Martin typeiS and a plane- wave basis-set 
with a cut-off energy of 30 Ry are used. The exchange- 
correlation potential is approximated within the LSDA in 
the parametrization of Perdew and Zunger^ In order to 
reduce the artificial defect-defect interaction, large super- 
cells with 216 lattice sites for 3C-SiC and 128 sites for 4H- 
SiC are used for all calculations of the defect energetics. 
For the 216-site cell, the Brillouin zone is sampled by the 
r point, whereas for the hexagonal 128-site cell a special 
Monkhorst and PacltiS 2x2x2 mesh is used. The values 
of the formation energies are given assuming silicon-rich 
conditions. Under carbon rich conditions they should be 
reduced by roughly 3Afff = 1.74 eV, where —AHf is the 
heat of formation of SiC. The vibrational properties are 
calculated using the frozen phonon method. The dynam- 
ical matrix 
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is calculated by applying displacements Qj from the equi- 
librium configuration and evaluating the force derivatives 
dFi/dQj. To obtain noticeable energy changes, the dis- 
placements Qj ought to be much larger than realistic 
phonon amplitudes. Therefore it is important to elimi- 
nate the anharmonic contribution when evaluating 
We do this by applying three different values of the dis- 
placement and extracting the linear force constant from a 
polynomial expansion of the forces. The forces have been 
converged to a relative accuracy of 10"'*. For 3C-SiC, 
the dynamical matrix is evaluated for a supercell with 64 
or 216 sites. The experimental bulk phonon modes are 
reproduced within 5 percent by this method. Yet the po- 
larization splitting of the TO and the LO modes at the 
r point is missing, since the macroscopic polarization of 
the crystal is incompatible with the periodic boundary 
conditions imposed on the supercell. We obtain a triple 
degenerate mode at 115.3 meV, in place of the TO mode 
of 98.7 meV and the LO mode of 120meV.i^ For 4H-SiC 
and for large clusters in 3C-SiC, the computational cost 
of evaluating the dynamical matrix of the whole 128-site 
or 216-site supercell is prohibitively large. In these cases 
we constrain the LVM calculation to the defect molecule 
(i.e. the dumbbell and its nearest neighbors) embedded 
in a supercell. We have examined the error resulting from 
the defect molecule approximation and from the uncer- 
tainty in the lattice constant. The results of these tests 
for a dicarbon antisite in 3C-SiC (see below) are shown in 
Fig. [3 The lines represent the frequencies of the LVMs 
versus the defect molecule size calculated at the theoreti- 
cal (LSDA) lattice constant. The frequencies shift within 
the grey area when the lattice constant is increased to- 
wards the experimental value of 8.239 Bohr. The black 
bars are the same frequency ranges, but in this case both 
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FIG. 2: Defect modes of (C2)si vs. the defect molecule size 
in a 216 sites cell. The straight line denotes the frequency 
at the theoretical lattice constant, the gray area is the fre- 
quency shift for the variation of the lattice constant between 
its theoretical and experimental value. The black bars are the 
corresponding values for the 64 sites cell as given in Tab. HI 
the arrows indicate the phonon replicas of the Dn centeri- 



the electronic structure calculation and the LVM calcu- 
lation have been performed in a 64 sites cell. We have 
also verified that under the pressure exerted by the peri- 
odic array of defects onto the lattice the lattice constant 
lies between the indicated limits. The LVM's variation 
due to the uncertainty in the lattice constant is thus less 
than 10 meV. It is also apparent that the defect molecule 
approximation affects the least localized modes 1 and (to 
a much lesser extend) 2. The largest effect on the accu- 
racy has the inclusion of the full supercell, which shifts 
the mode 1 by 14meV in the 216 sites cell. Compared to 
the full 64 sites cell, the mode 1 is only 2 meV higher. 



III. THE DICARBON ANTISITE 

Let us now consider one by one the candidates for 
the Dn center as listed above. First, we find that the 
originally suggested carbon di-interstitial has a forma- 
tion energy of about 12 eV. This is extremely high even 
for SiC and practically rules out this model. A simi- 
lar defect is the carbon split-interstitial Csp{ioo)j i-e. a 
pair of carbon atoms sharing the same (carbon) site and 
pointing in (100) direction (in 3C-SiC). This defect also 
possesses a short carbon-carbon bond and it is energet- 
ically the most favorable carbon interstitial.^"' The sec- 
ond lowest carbon interstitial is the carbon-silicon split- 
interstitial Cgpgi^ioo) (a carbon and a silicon atom on a 
silicon site) . Here the additional carbon atom is also con- 
nected to its nearest neighbors by short carbon-carbon 
bonds. However, according to our calculations both de- 
fects are very mobileyi^ which is in conflict with the ther- 
mal stability of the Dn center. In addition, we find the 
LVMs of both defects are not compatible with the Dn 
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spectrum. The neutral carbon split-interstitial Csp{ioo) 
in the tilted structurei^ has two LVMs at 120.4 and 
188.7 meV and LVMs in the phonon gap. The carbon- 
silicon split-interstitial C^^gj^j^pg^ possesses a richer vibra- 
tional spectrum but still with only three non-degenerate 
LVMs above the bulk spectrum. This also does not match 
the number of the characteristic Dn lines. Another possi- 
ble candidate could be the carbon antisite Cgi- Following 
the carbon vacancy Vc it is the second most abundant 
defect in SiC. However, it shows only vibrational reso- 
nances in the bulk phonon spectrum and has no true 
LVMs. Also small clusters of carbon interstitials do not 
provide LVMs compatible with the Dn spectrum.^ 

The simplest candidate that possesses the main fea- 
tures of the Dn center is the dicarbon antisite (C2)si 
that we have already mentioned above in connection with 
Fig. [3 The dumbbell of two carbon atoms sharing a sil- 
icon site (highlighted atoms in Fig. ^ left) is surrounded 
by a tetrahedron of four carbon atoms. In 3C-SiC the 
dumbbell is pointing in (100) direction whereas in 4H- 
SiC it is oriented either along the (1011) or the (1011) 
direction. Depending on the Fermi level, the defect may 
acquire the charge states from 2+ through 2^. 

We find that the dicarbon antisite has a relatively low 
formation energy. For the neutral charge state the val- 
ues are 6.7 eV and 7.8 eV in 3C and 4H polytypes, re- 
spectively. In 4II-SiC two possible substitutional config- 
urations exist, i.e. hexagonal (h) and cubic (k), but the 
formation energy at these two sites only differs by 0.1 eV. 
These values are comparable to the formation energies of 
Csp{ioo> and CspSi(ioo>- Still, the calculated formation 
energy indicates that its equilibrium abundance is too 
low for experimental observation. Hence kinetic effects 
should play an important role in the dicarbon antisite 
formation. On the other hand, once created, the defect 
has a very high stability, as reflected by the high bind- 
ing energy. Before we turn to kinetic aspects, we discuss 
the electronic properties and the LVMs of the dicarbon 
antisite. 

In 3C-SiC, the dicarbon antisite possesses two degen- 
erate defect orbitals with energy levels in the band gap. 
In the neutral charge state these orbitals are occupied 
by two electrons, which can form either a spin (low- 
spin) or a spin 1 (high-spin) configuration. According 
to the calculation, the high-spin state is about SOmcV 
lower in energy than the low-spin state. This energy dif- 
ference is, however, too small to uniquely determine the 
ground state spin configuration. For this reason we con- 
sider both options. The defect wave function consists of 
bonds between the three upper carbon atoms and the 
three lower carbon atoms of the defect tetrahedron (cf. 
Fig. n left). In the low-spin configuration, one of these 
bonds is fully occupied, resulting in a splitting of the 
defect levels and a Jahn- Teller distortion, with the up- 
per or the lower pair of the tetrahedron being stretched. 
This reduces the original D2d symmetry to C2v In the 
high-spin configuration, both bonds are equally occupied, 
retaining the degeneracy and the D2d symmetry. As we 



discuss below, this geometrical difference is reflected in 
the vibrational modes of the defect. In 4H-SiC, the high- 
spin state is about 120 meV lower than the low-spin state, 
which is — as in 3C-SiC — a very small energy difference. 
Due to the lower symmetry of the crystal, the symmetry 
group is always Cih- The defect's structure on the cubic 
and the hexagonal site is qualitatively similar. 

As outlined above, the dicarbon antisite is the simplest 
defect that possesses LVMs that are similar to the phonon 
spectrum of the Dn center. The calculated frequencies 
(using a 64-sites supercell) for the neutral charge state in 
low-spin and high-spin configurations are hsted in Tab.|l| 
along with the experimental data for the five character- 
istic Dn phonon modes and our calculated results for the 
two neighboring dicarbon antisites. We see that in the 
3C polytype the low spin configuration provides an LVM 
pattern similar to the five phonon replicas of Dn- Similar 
results are obtained for 4H-SiC, except that the frequen- 
cies of the modes 1 to 4 are lowered due to the defect 
molecule approximation. 

In spite of a different crystal symmetry, the vibrational 
patterns of the LVMs in 3C- and 4H-SiC are very similar, 
owing to the localized nature of the vibrational modes. 
The highest mode 5 is a stretching vibration, where the 
two atoms of the dumbbell oscillate against each other. 
Its calculated frequency is about 8% higher than the ex- 
perimental value. However, this is the usual accuracy of 
these defect calculations.^^ The modes 3 and 4 represent 
an oscillation of one of the dumbbell atoms against its 
neighbors. In 3C-SiC, these modes are affected by the 
Jahn- Teller distortion of the low-spin configuration. As 
expected from the symmetry, the modes are degenerate 
in the high-spin and non-degenerate in the low-spin case. 
This effect can also be seen in Fig. |2 where the 216 
sites calculation, although performed spin-unpolarized, 
has been performed in a high-spin like configuration. 
Considering the small energy difference between the high- 
spin and the low-spin configuration both options should 
be allowed. In addition to the already mentioned clus- 
tering of carbon atoms, this may be a source of further 
lines in the luminescence spectra. In 4H-SiC, the modes 
3 and 4 are always non-degenerate, although the split- 
ting in the high-spin configuration is smaller than in the 
low-spin configuration, due to the lower symmetry of the 
crystal field. The mode 2 describes the vibration of the 
dumbbell along its axis against the enclosing tetrahedron. 
The mode 1 is the least localized mode and therefore it is 
most affected by the employed supercell technique. This 
mode corresponds to a breathing-like vibration of the 
enclosing tetrahedron. While for the low-spin case the 
breathing mode lies above the highest bulk vibration, it 
is in resonance with the bulk spectrum and lies 8 meV 
below the calculated value of the highest bulk mode for 
the high-spin case in a 64 sites cell. This is in agreement 
with Gah et al. (Ref. El) who considered the high-spin 
state of (C2)si- In 4H-SiC, a similar behavior is not obvi- 
ous from our results, as the effect may be masked by the 
defect molecule approximation. In addition to the high- 
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TABLE I: LVMs of (C2)si and ((02)31)2 in the neutral charge state in meV calculated at the LSDA lattice constant. Is denotes 
the low-spin, hs the high-spin state. The subscripts h and k refer to the hexagonal and the cubic sites. For 4H-SiC and 
((02)31)2 the defect molecule approximation has been used. For the Dn center only the five highlighted frequencies are given 
(cf. Refs. 2 7). 



LVM 3C 4H 





(C2)si, Is 


(C2)si, hs 


((C2)si)2 


exp. 


(C2)3i,h, Is 


(C2)si,h,hs 


(C2)si,k,ls 


(C2)3i,k,hs 


((02)31)2, hk 


exp. 


1 


116.4 


107.6 


123.3 


127.8 


103.0 


100.2 


102.3 


101.5 


119.6 


127.1 


2 


128.5 


123.2 


125.6 


130.9 


121.9 


121.0 


119.7 


120.7 


132.2 


129.8 


3 


139.4 


141.7 


127.1 


146.3 


136.2 


136.4 


135.0 


136.2 


147.0 


146.1 


4 


148.3 


141.8 


139.3 


152.2 


139.3 


137.6 


139.1 


137.6 


160.2 


152.4 


5 


177.5 


177.5 


168.6 


164.7 


179.4 


177.1 


178.0 


177.0 


162.4 


164.4 


6 






169.7 

















frequency LVMs discussed above, we also find maxima of 
the vibrational density of states in the phonon band gap. 

It is noteworthy that the resuhs do not change sig- 
nificantly for different possible charge states of the 
deiectm^^ We observe a slight lowering of the LVM's fre- 
quencies with the increase of the number of localized elec- 
trons, resulting from an expansion of the defect molecule 
due to the occupation of the bonds. For 3C-SiC, the 
modes 3 and 4 are degenerate in the charge states 2+ and 
2^ , since the electronic defect levels are either empty or 
fully occupied, and a Jahn- Teller distortion is suppressed. 

IV. CLUSTERING OF CARBON ATOMS 

Usually, the Dh center is observed in ion-implanted 
and annealed samples. It is stable up to temperatures 
of 1700°C. This and the previously discussed low equilib- 
rium concentration suggest a kinetically controlled for- 
mation of this defect from non-equilibrium vacancies and 
interstitials. At first, a carbon antisite is formed via 
a recombination of a carbon interstitial with a silicon 
vacancy. It merges then with a mobile carbon split- 
interstitial, producing a dicarbon antisite with a high 
binding energy that ranges from 3.9eV ((C2)g^) to 5.1 eV 
((C2)g^). The dicarbon antisite can serve as a con- 
densation center for larger clusters. Such clusters grow 
by aggregating carbon split-interstitials around the core. 
Adding a further carbon atom to the dicarbon antisite 
releases a binding energy of about 4.8 eV. Binding one 
more carbon interstitial to this complex yields a very 
small cluster of four carbon atoms at a silicon site (cf. 
Fig- n right). According to our calculations the energy 
of 2.8 eV is needed to remove a carbon atom from this 
structure. However, for this defect we obtain LVMs well 
above the Dn spectrum. A detailed discussion of carbon 
clusters will be given elsewherCiiS; 

In addition to the clusters located around a single 
silicon site, aggregates of several dicarbon antisites are 
also possible. A pair of the dicarbon antisites in 3C- 
SiC ((€2)31)2 is shown in Fig. In 4H-SiC three dif- 
ferent configurations, which can be compared to its 3C- 




FIG. 3: Structure of a pair of dicarbon antisites ((02)31)2 in 
30-SiO. 

SiC counterpart, are possible when the two di-antisites 
occupy two cubic (fcfc), two hexagonal (hh) or cubic- 
hexagonal (fc/i) sites. As an exemplary structure we have 
considered in detail the fc/i-configuration, which has a 
local geometry similar to the 3C-SiC case (Fig. but 
the bond length between the neighboring dicarbon anti- 
sites is shorter. We found that this defect is very stable 
in both polytypes: in 3C-SiC it costs about 5.9 eV to 
remove a carbon atom, whereas in 4II-SiC this energy 
is 6.7 eV. The calculated LVMs of both complexes are 
given in Tab.^ They have been calculated using the de- 
fect molecule approximation that allows the vibration of 
the two dumbbells and of the surrounding carbon neigh- 
bors within the 216 (3C-SiC) and 128 (4H-SiC) sites ceU. 
As discussed above, the numerical results are affected by 
the defect molecule approximation. Alternatively, one 
could work with the smaller 64 sites cell which would al- 
low to calculate the dynamical matrix of the whole cell. 
However, we find that in this cell the electronic structure 
is so strongly distorted by the proximity effects that it 
is unacceptable for the analysis of the LVMs. For the 
pair of dicarbon antisites the highest almost degenerate 
modes 5 and 6 involve a stretching vibration of the dumb- 
bells, similar to the lone dicarbon antisite. Due to the 
interaction of the dumbbells these modes become softer 
and combine into anti-symmetric and symmetric vibra- 
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tions. The modes 3 and 4 represent a motion of the 
two lower and upper carbon atoms of the dumbbeUs in 
bond direction against their surrounding neighbors, re- 
spectively. In 4H-SiC, the defect shows a similar be- 
havior, but the energy of all modes is by about 8meV 
lower than in 3C-SiC so that the lowest mode (with an 
energy of 114.4 meV) drops into the bulk phonon spec- 
trum. This effect and also the larger frequency differ- 
ence between the stretching modes reflects the stronger 
relaxation of the defect in this polytype. We have re- 
cently learned that A. Gali et al. calculated the LVMs 
of the /c/c-configuration, finding results similar to the 3C- 
SiC values, and speculated that this defect was the Dn 
centeriSi As the fcfc-configuration has a 3C-like environ- 
ment (in contrast to the /c/i-configuration), the similar- 
ity to 3C-SiC is quite natural. However, the conjecture 
that the ((02)31)2 model completely describes the Dn 
center does not seem to be supported by experimental 
data. Recent experiments have investigated the multi- 
tude of zero phonon lines (ZPLs)^ as well as the phonon 
replicas^ of the Dn center. In 6H-SiC, Sridhara et al^ 
concluded from temperature and stress-dependence ex- 
periments that the four ZPLs of Dn stem from the dif- 
ferent excitonic states of the bound exciton (one ground 
state and three excited states), but not from different 
binding centers of the same exciton. Regarding the 
phonon replicas in 4H-SiC more than 12 LVMs and a 
single ZPL have been reported by Carlsson et al.,^ which 
is consistent with the earlier results for 4H and 6H-SiC 
of Sridhara et alX In addition, as we noted above, a com- 
parison of the spectra in 411 and 6H-SiG- indicates 
that more than five LVMs may be present also in 3C-SiC. 
The apparent large number of phonon replicas in 411 and 
6H-SiC makes it tempting to assume that they originate 
from different geometrical configurations of the same sim- 
ple defect. For example, for the ((02)31)2 in 4H-SiC this 
would be kk, hh and kh. This assumption has, how- 
ever, the following implications. Firstly, it follows that 
a larger number of LVMs should be present in 6H-SiC 
than in 4H or 3C-SiC, since less or no inequivalent sites 
are available in the two latter polytypes. Secondly, the 
different configurations should give rise to a multitude 
of exciton ground-state ZPLs. According to the conclu- 
sion of Sridhara et al^ and the rich vibrational structure 
contained in all polytypes (cf. discussion above) this is 



not evident from the experiments. This applies, in par- 
ticular, to the ((€2)31)2 defect. Furthermore, the LVM 
pattern of this defect is, in fact, not compatible with 
the Dn spectrum. The small splitting of the two highest 
LVMs (see Tab. Pi, which persists through different poly- 
types and configurational models, is not observed in the 
Dn spectrum which has nearly equidistant lines. Thus 
the ((€2)31)2 model alone cannot explain all features of 
the Dn center. Larger carbon clusters are likely to be 
involved. For such clusters the differences between the 
polytypes lose their importance. This may be responsi- 
ble for the polytype-independent LVMs of the Dn cen- 
ters. Further experimental data are necessary to clarify 
the situation, especially for the 3C polytype, which lacks 
inequivalent lattice sites. 



V. CONCLUSION 

In conclusion, we have analyzed the structural, elec- 
tronic and vibrational properties of the dicarbon antisite 
complex. The defect possesses a characteristic vibra- 
tional pattern in the frequency range of the Dn spec- 
trum, but cannot explain all the phonon replicas ob- 
served. Yet, we find that the dicarbon antisite can act as 
an aggregation center for carbon clusters that allow for a 
much richer phonon spectrum. All the defect complexes 
presented above possess very high dissociation energies, 
which guarantee their high thermal stability. Due to the 
high formation energies, these defects can exist in as- 
grown material only in very low concentrations but can 
be kinetically created during the annealing. Yet, none of 
the considered defect complexes can explain all features 
of the Dn center. Most probably the Dn-type centers 
are related to larger carbon clusters with the dicarbon 
antisite playing the role of an elementary building block 
of the aggregates. 
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